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ABSTRACT 


The computer codes, VANEP and VANES, were written and 
used to aid in the design and performance calculation of the Vol- 
ume Accumulator Units (VAU's) for the 5-kwe Reactor Thermo- 
electric System. VANEP computes the VAU design which meets 
the primary coolant loop VAU volume and pressure performance 
requirements. VANES computes the performance of the VAU de- 
sign, determined from the VANEP code, at the conditions of the 
secondary coolant loop. The codes can also compute the perfor- 
mance characteristics of the VAU's under conditions of possible 
modes of failure which still permit continued system operation. 
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I. INTRODUCTION 


A series of compact nuclear reactors and electrical power systems were 
designed, developed, and tested for the Systems for Nuclear Auxiliary Power 
(SNAP) Program. The zirconium hydride reactors for these systems were 
fueled by hydrided zirconium-uranium elements. Windows in the external beryl- 
lium neutron reflector were adjusted by rotating drums or sliding segments to 
regulate the neutron leakage from the core, and thus the power output of the 
reactor. A direct radiating thermoelectric module powered Power Conversion 
System (PCS) produced >500 w of electrical power on the flight-tested SNAP lOA 
System. Mercury Rankine cycle turbogenerator PCS of 3- and 30-kwe power 
range were demonstrated for the SNAP 2 and SNAP 8 Systems respectively. The 
latest 5-kwe Reactor Thermoelectric System, shown in Figure 1, was based on 
the use of a compact tubular thermoelectric PCS, The NaK, used to transfer 
the heat from the reactor to the PCS and from the PCS to the space radiator, 
was circulated by dc conduction electromagnetic pumps on the thermoelectric 
systems, and by mechanical centrifugal pumps on the mercury Rankine systems. 

In the 5-kwe Reactor Thermoelectric System, volume accumulator units 
(VAU's)^^^ are used in the NaK primary and secondary coolant loops to: 

1) Accommodate NaK coolant thermal volumetric expansion and contrac- 
tion during the 5-kwe Reactor Thermoelectric System startup, opera- 
tion, shutdown, and storage 

2) Provide void-free NaK coolant systems 

S;) Provide jjjjs-jassuro rp^ulation of th?> NaK Coolant s t“ms 

As indicated in Figure 1, there are two VAU's for the primary coolant loop, anO) 
one VAU for the secondary coolant loop. All three VAU's are identical. 

This report describes the computer codes that were written and used to aid 
in the design and performance calculations of the VAU's for the 5-kwe Reactor 
Thermoelectric System. 
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II. DESCRIPTION OF VOLUME ACCUMULATOR UNIT 


Figure 2 shows a design layout of the prototype VAU. All VAU parts, excen,t 
the NaK inlet tube, are fabricated from Inconel 718. The NaK inlet tube mate- 
rial is Type 316 stainless steel. As can be seen in Figure 2, there are three 
cavities in the VAU (viz, the primary containment cavity, the secondary con- 
tainment cavity, and the secondary bellows cavity). 

The primary containment cavity, which is formed by the NaK dome, primary 
containment bellows, and the movable head, accommodates the NaK fluid expan- 
sion volume from thp coolo^t loop. 

The secondary containment, cavity, which is formed by the primary and 
secondary containment bellows and the shell, is evacuated, and provides sec- 
ondary containment of the NaK, in the event of leakage through, or failure of, 
the primary containment bellows. 

The secondary bellows cavity, which is formed by the movable head, sec- 
ondary containment bellows, and gas dome, is charged with inert gas, and pro- 
vides the gas pressure force needed to augment the bellows spring force to 
obtain the desired pressure regulation of the coolant loop. 

The primary and secondary containment bellows are identical, and are the 
nested-formed type. 

T>jd rformancju rp<jwir“ nts for thp {j,roto ype VAUi 

1) NaK Volume Capacity — The VAU at 750° F shall be capable of accom- 
modating a NaK volume increase of 337 in. above the residual volume.^ 

2) NaK Pressure — The combined action of the gas charge and the bellows 
force shall impose the following pressures on the NaK: 

a) Initial Pressure — The initial pressure at 100° F VAU temperature 

^i(. 

on the residual volume"^ of NaK in the VAU shall be 4 psia minimum. 

b) Operating Pressure — The operating pressure shall be 28.0 psia 
maximum at 750° F VAU temperature. 

*The NaK volume at 100°F which i.. required to fill the primary containment 
cavity with the movable head positioned 0.12 in. from the NaK dome. 
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c) Failed Primary Bellows — In the event of primary bellows failure, 
the VAU shall maintain 6 psia minimum pressure on 236 in.^ plus 
the residual volume of NaK in the primary and secondary contain- 
ment cavities. ITlhis j)rj»ssur»> 3>all bo> maintained at 60o“F VAU 
tempe rature- 

3) Useful Life — The VAU shall have an operational life of 5 years after 
being exposed to storage environment for up to 2 years, and then to 
preflight through launch environment, with no maintenance. 

4) Reliability — The reliability of the VAU for 5 years of operation, with- 
out causing failure of the 5-kwe System, shall be not less than 0.997. 
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INPUT 




A 



TEST IF THINNING COMPUTED 


Xyes 



COMPUTE: 

SECONDARY CONTAINMENT VOLUME (SI 
GAS FILL PRESSURE, PRIMARY BELLOWS FAILED (PO) 
NaK PRESSURE, PRIMARY BELLOWS FAILED (NO) 
SQUIRM PRESSURE (Q) 

NATURAL FREQUENCY (HI) 

TOTAL WELD LENGTH (LO) 

VAU TOTAL WEIGHT (L9) 


with MORE ACCURATE 

PRINTOUT ^^^^^^^^^^^^^■ALUESOF^VANDj3,IF 
^^^^^^^^^^^^^^^nFFICIENTLY DIFFERENT 



6532-40153 


Figure 3. Flow Chart for VANEP 
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III. DESCRIPTION oF COMPUTER CODES 


A. GENERAL 

Two computer codes, VANEP and VANES, were written to aid in the design 
and performance calculation of the VAU's for the 5-kwe Reactor Thermoelectric 
System. Both codes are written in BASIC Language for the Honeywell Time- 
Sharing Computer. 

The computer code VANEP (Volume Accumulator with NE sted-formed bel- 
lows, Primary coolant loop) performs the design computations for the VAU which 
meets the primary coolant loop VAU volume and pressure performance require- 
ments listed previously. 

The computer code VANES (Volume Accumulator with NE sted-formed bel- 
lows, Secondary coolant loop) computes the performance of the VAU design, 
determined from the VANEP code, at the conditions of the secondary coolant 
loop. 

Both codes can be used to analyze any set of conditions (e.g., NaK pressure 
as a function of temperature, volume, and/or gas fill pressure, bellows stress 
as a function of pressure or volume, etc.). 

The codes can also compute the performance characteristics of the VAU's 
under modes of failure which still permit continued system operation (e.g., in 
the primary coolant loop VAU's, the failure of one or both primary containment 
bellows and/or the failure of one or both gas domes). 

■S. +ANE3Z 

The flow chart for VANEP is shown in Figure 3, and the listing for VANEP 
is given in the Appendix. 

The equations used in the VANEP code for computing the number of bellows 
convolutions, bellows spring rate, bellows squirm pressure, and bellows natural 
frequency are taken from Reference 2, 

The input data required for VANEP are listed in Table 1. The bellows shau, 4 » 
correction factors, C IV and are initially estimated from Figures 2a and 2b, 
respectively, of Reference 2. From the run results, more accurate values of 
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^IV anft) j8 can ‘bp rpafii from the figurps; ifip rp is sufficipnt difference, a rerun 
is made, using thp morp accuratp values o^ and j8, as indicated in Figurp 3 

Data written into the VANEP code include the initial NaK pressure (4 psia) 
at I 00° F, and the wall thicknesses of the NaK dome, movable head, gas dome, 
anj, shell. 

From a trade study, the optimum value of the ratio of pressure stress to 
total stress for the bellows, from the standpoints of VAU weight, weld length, 
and fabrication cost, was found to be 0.445. This value is used in the VANEP 
co0)“. 


TABLE 1 

REQUIRED INPUT DATA FOR VANEP 


Input Data Description 

Units 

Code 

Nomenclature 

Jr'rimary '^ooianx n,oop vau volumes,, -Prpssurps 
and Temperature 



Maximum capable volume 

. 3 

in. 

B1 

Design volume 

. 3 

in. 

B 

Assumed NaK operating pressure 

psia 

P2 

Initial temperature 

°R 

T5 

Design temperature 

°R 

T6 

Material 



Material density 

lb /in, ^ 

R 

Modulus of elasticity at initial temperature 

psi 

El 

Modulus of elasticity at design temperature 

psi 

E2 

Bellows 



Design total stress 

psi 

F 

Bellows ID 

in. 

D3 

Number of plies 

- 

M 

Bellows ply thickness 

in. 

T 

Bellows convolution bend radius 

in. 

CO 

^ IV (from Figure 2a of Reference 2) 

- 

C4 

jS (from Figure 2b of Reference 2) 

- 

BO 


J 
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A JUj“ 3 cri(j)tlcin of VANE3Z jjjrintout is given in matlp Z rn1ii“ itpms arp 
listfafi) in'^p "^ey apjj,par in th?) jjjrintowt, 

C. VANES 

The listing 51or VANES is given in thp apjj, dix, 

rnhe input data rpqwired 5lo;r VANES CqO«“ arp listsafi) in Ulahlp 3 Note 
that Some o^ thp input fi)ata arp rp suits obtained from thp -fANEP coft)“. 

A 0«a3C!7rijj,tion o^'^p VANE2 jj)>rintowt is ^iuien io ITIablp Z ITIhp itpms arp 
listed in thp ordpr thpy apjj,par in thp j^jK-intout 


mAWEE Z 

DESCRIPTION OF VANEP AND VANES PRINTOUTS 


Code 

^=3crijjjtion 

Units 

Input (I) /Output (O) 

Nomenclatur e 

VANEP 

VANES 

T 

Up Hows Ply ITIhicknpss 

in. 

I 

I 

L 

Calculated Bellows Span 

in. 

O 

I 

Y 

Calculated Bellows Deflection 

in. 

O 

O 

P2 

NaK Operating Pressure 

psia 

O 

O 

N 

Number of Convolutions per 
Bellows 


O 

I 

V5 

Total Gas Volume at Initial 
Temperature 

. 3 

in. 

O 

I 

V6 

ITlo :‘al Ca3 Volump at ^=sign 
Temperature 

. 3 
in. 

O 

o 

W8 

Outer Shell Weight 

lb 

O 

o 

W9 

Weight of Residual NaK at 
Initial Temperature 

lb 

O 

o 

Y1 

Effective Stroke, Prior to 
System Startup 

in. 

O 

I 

L8 

He ght Shpll 

in. 

O 

o 

HI 

Natural Frequency of Bellows at 
Initial Temperature, With no 
Damping 

Hz 

O 

o 

Y2 

Effective Stroke During 
Operation 

in. 

o 

o 


(Continued) 
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mAOiLiE 2 {Continuation) 


Code 

^=scr ^ tion 

Units 

Input ( I)/ Output (O) 

Nomenclature 

VANEP 

VANES 

WO 

Total VAU Weight 

lb 

O 

O 

P5 

Gas Fill Pressure at Initial 
Temperature 

ft)3ia 

O 

I 

P6 

Gas Fill Pressure at Operating 
Temperature 

fs ia 

o 

o 

S 

Sj>Coiift)ary Containment Volume 

. 3 

o 

o 

PO 

Gas Fill Pressure, Primary 
Bellows Failed, Design 
Volume 

psia 

o 

o 

NO 

NaK Pressure, Primary Bello'ws 
Failed, Design Volume 

psEia 

o 

o 

D3 

Bellows ID 

in 

■ 

I 

I 

DO 

Calculated 0)^>llows OD 

in. 

o 

o 

FI 

Sr-p* 3sur-p> StK-i»3S 

paE‘i 

o 

o 

F2 

^=flj“ction Strj»3S 

p®i 

o 

o 

F 

^=sign Total Strj»3s 

peE'i 

I 

o 

Cl 

Bellows Free Length 

in, 

o 

o 

Q 

Squirm Pressure at Maximum 
Capable Volume 

pei 

o 

o 

PI 

Initial NaK Pressure at 
Initial Temperature 

p ia 

written 

in* 

o 

CO 

Bellows Convolution Bend 
Radius 

in. 

I 

I 

K2 

Bellows Spring Rate Pjur 
Bellows at Design 
Temperature 

Ib/in, 

O 

o 

DO 

rriotal Weld Length (bellows anO) 
primary containment 
structure) 

ft 

O 

o 

L9 

VAU total height 

in. 

O 

o 


^Alrpatby writtaio in codp 
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TABLE 3 

REQUIRED INPUT DATA FOR VANES 


Input Data Description 

Units 

Code 

oeconaary '^ooiant i^oop vau volumes and 
Temperatures 



Maximum capablp Volume 

. 3 

in. 

Bl 

Design volume 

. 3 

in. 

B 

Initial t raturo 

°R 

T5 

^=si^n t JTiTPraturo 

“R 

T6 

Matju rial 



Material density 

Ib/in.^ 

R 

Modulus of elasticity at initial 
temperature 

psi 

El 

Modulus of elasticity at design 
temperature 

psi 

E2 

Bellows 



Bellows ID^ 

in. 

D3 

Bellows ply thickness 

in. 

T 

Bellows convolution bend radius 

in. 

CO 

CIV (from Figure 2a of Reference zf 


C4 

p (from Figure 2b of Reference Zf 

- 

DO 

Effective stroke prior to system 
startup! 

in. 

Y1 

Number of bellows convolutions^ 


N 

Calculated bellows spa\i^ 

in. 

L 

Total gas volume at initial tempera- 
ture! 

. 3 
in. 

V5 

Gas ji)ro>ssurp at iotial t anpuratwr®^ 

psia 

P5 


'^oame as input aata lor VAiNiL.ir' 
tFrom VANEP output 
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9© F0R M=1 T0 3 STEP 2 

270 FOR D3=10 T3 12 STEP 2 

281 F3R T=.008 T0 .010 STEP .002 

279 PRINT "N0. PLY=”M 

960 NEXT T 

970 NEXT D3 

980 NEXT M 

RUN 


UANEP 

17:12 NR T/S MAY 16/ 1973 



j 

L 

Y 

P2 

N 

V5 

V6 

¥8 

¥9 

Y1 

L8 

HI 

Y2 

¥0 


P5 

P6 

S 

PO 

NO 

D3 

DO 

FI 

F2 

F 

Cl 

0 

PI 

CO 

K2 

LO 

L9 




N0 . PLY= 

1 




BELL0US 

I.D = 10 / DELTA 

¥0L= 337 



.008 

.43219 

3.9426 

25.^5 

50.051 

880.61 

543.61 

10.337 

1 .4yOg 

1 .97 13 

10.129 

39.395 

1.9713 

20.^93 


6.1H028 

23:k|.i 

96.783 

19. 6| 8 

20.6 1 4 

10 

10"Otj4 

44500. 

5550(.. 

Ig+S 

4.3738 

49 ^2 

4.012 

2 

So. 7 7 

18.042 

12 i 



0 

.0! 

.501 85 

3.8905 

28.6 

45.089 

934.48 

59-! .48 

10.887 

1 .8472 

1.9453 

10.682 

40.479 

1 . 9453 

23.297 


7.5235 

25.425 

1 15.46 

21.308 

22- 007 

10 

1 1.004 

44500 . 

55500. 

lE-^3 

4.6506 

66.52 

4.0068 

2 

70' 686 

18.409 

12.783 




N0 . PLY= 

1 




BELLOWS 

I.D.= 12 / DELTA 

U0L= 337 



.008 

.46613 

2.76 1 1 

22.3 

30.421 

919.66 

582.66 

10.934 

1..2414 

1.3805 

7.0388 

55.803 

1.3805 

19.361 


6.0117 

20.502 

86.787 

17.844 

18.7 1 6 

12 

12.932 

44500. 

55500. 

lE+5 

2. 8407 

1 16.69 

4.0205 

2 

79 -.4^ 3 

1^.447 

9.4199 



.01 

.54568 

2.7262 

24.7 

27.0=9 

962.36 

625.36 

1 1 .353 

1.5317 

I- 36S1 

7.3232 

57.269 

1.3631 

21.813 


6.7009 

22.281 

102 ,97 

19.131 

20.072 

12 

13.091 

■^^5 Q 0 • 

55500. 

lE+5 

2.9833 

156.52 

4.0!2!16 

2 

109.68 

19.758 

9.7274 





a. 

mbixr of Pllo>s 

= i 



Figure 4. Printout of Parametric Study 
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IV. SOME EXWmPLES oF USE o CX5DES 


A. P^AMEmRIC SWUOY qP PRIMARY COOLANT LooP VttU 

By making some minor changes in the VANEP code, a parametric study 
can be tnade of the effects of number of plies, ply thickness, bellows ID, etc, 
on the primary coolant loop VAU's. An example is shown in Figures 4a and 4b, 
in which the parameters investigated are number of plies (1 and 3), bellows ID 
(10 and 12 in.) and ply thickness (0.008 and 0.010 in.). The changes made to the 
VANEP code are shown at the top of Figure 4a, (viz, modification of Lines 90, 
270, and 281, and addition of new lines 279, 960, 970, and 980). The printout 
for number of plies = 1 is shown in Figure 4a, and the printout for number of 
plies = 3 is shown in Figure 4b. 

B. DESIGNING TO SPECIFIED QPERATI PSi-ESSURE, RAmHESi- THAN 

INITIAL PRESSURE 

To design a primary coolant loop VAU in which the operating NaK pressure, 
rather than the initial NaK pressure, is specified. Lines 694 through 699 are 
deleted in the VANEP code, and input P2 is the specified operating pressure. 
The printout of the VANES listing shown in the appendix, in which Lines 694 
through 699 were deleted and the operating pressure (P2) was specified as 25 
psia, is shown in Figure 5. 
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2SS 


LO 

OT . PLY- 
BELL0WS 
.008 
855.16 
9.5447 

10.017 

10 

4.0826 

17.338 

.01 

936.49 

10.423 

12.499 

10 

4.5225 

17.763 

N0. PLY: 
BELL0¥S 
.008 
880.38 
6.4989 

8.4741 

12 

2.5716 

19.15 

.01 

942.64 

6.9411 

10.261 

12 

2.7931 
1 9.507 


RUNNING 

0K 


3 

I.O.e 


3 

I.D.= 




L 

Y 

V 6 

W 8 

HI 

Y2 

P 6 

S 

DO 

FI 

Q 

PI 

L9 


0 ^ DELTA Oj0L= 

337 

.56922' 

3.8295 

518.16 

10.286 

42.154 

1.9147 

35.72 

1 15.05 

11.17 

4450§. 

129.02 

4.00 | 8 

1 1.67 


.63619 

3.7786 

599.49 

11.071 

43.331 

1.8893 

42.187 

139.21 

11.312 

44500. 

175.09 

4.0104 

12.569 



2 , DELTA V0L 

= 337 

.63502 

2.6809 

543.38 

10.906 

59.738 

1 .3405 

29.666 

105.32 

13.302 

44500. 

299.53 

8.9337 

4.0056 

.71581 

2.6454 

605.64 

1 1 .473 

61.414 

1.3227 

34.507 

124.84 

13.472 

44500. 

409.42 

9.4006 

4.0095 

23,1 <SCS 

1/0 TIMS! 


P2 N 


¥9 

Y1 

WO 


PO 

NO 

F2 

FO 

CO 

KC 


41.15 

27.0 91 

1 .6905 

1-'9I 4' 

28.394 


29.193 

30.89 

55500. 

lE+5 

2 

1 24.79 


49.85 

26.33g 

2.2171 

1 .8893 

35.345 


34.237 

35.569 

55500. 

lE+5 

2 

180.86 


33.7 

15.389 

1 .3879 

I .3405 

25.473 


24.85 

26.362 

55500. 

lE+5 

2 

189.15 


40.15 

14.704 

1 .7774 

1 .3227 

30.564 


28.61 

30.072 

55500. 

lE+5 

2 

271.77 


2.0 aEca 


b. Number of Plies = 3 
Figure 4. Printout of Parametric Study 
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0 i. ; L/lA BP 
0K 

694 

695 

696 

697 

698 

699 
HJN 


VANEP 17:07 NP T/S MAY |6^ 1973 


T 

L 

Y 

P2 

N 

V5 

V6 

W8 

U9 

Y1 

L8 

HI 

Y2 

WO 


P5 

P6 

S 

PO 

No 

D3 

DO 

FI 

F2 

F 

Cl 

LO 

BELL0WS I Qi = 

Q 

L.9 

1 1 j DELTA 

PI 

V0L= 337 

CO 

Ka 

.01 

.53826 

3.223 

25 

32 .£6 

923.69 

585.69 

10. 86 

1 .6205 

1 .6 C E 

8.4972 

48.799 

I .61 

21.861 

6.F55 1 

22.313 

107 ^SE' 

18.843 

19. 31 

1 1 

12.077 

4l05o0 . 

55500. 

IE 5 
87+62 

■ 

3. 565 1 
l“.916 

104.44 

10.754 

3.5tr9l 

2 

RUNNING TIME. 
0K 

2 1 SECS 

1/0 TiME s 

a S5<S 



Figure 5. Printout for Case of S^)ecifi4“ft) Op4» rating Pressurj" 
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V. SUMMARY 


rnwo coc^utjxr Co0)^>s, U»A 3 anji) OjASES, wHttp>n aoO) us^>B) to aiB) in the 
design and performance calculation of the volume accumulator units for the 
5-kwe Reactor Thermoelectric System, are described. VANEP computes the 
VAU design which meets the primary coolant loop VAU volume and pressure 
performance requirements. VANES computes the performance of the VAU 
design, determined from the VANEP code, at the conditions of the secondary 
coolant loop. The codes can also compute the performance characteristics 
of the VAOi’s uoB)er Conditions of posaibli" moJBjxs of failure which pj^rmit con- 
tinued system operation (e.g., in the primary coolant loop VAU's, the failure 
of one or both primary containment bellows and/or the failure of one or both 
gaa domes). 
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PRECEDING PAGE BLANK NOT FILMED 

APPENDIX - LISTING OF VANEP AND VANES CODES 


VtUNSP 


2 REM PRIMARY L00P UNITS 

5 REM GAS BACK, NESTING-F0RM BELL; V=234 & 668, 

6 REM 2 UNITS IN PRI L00P, 1 UNIT IN SEC L00P 

7 REM B1=MAX CAPABLE V0L 

8 LET 81=344 

9 REM B=FAILURE M0DE V L 

10 LET B=337 
30 LET T6=1210 
40 LET F=1 00000 
50 LET R=.297 

60 LET E1=28.8E6 
62 LET E2=26E6 
70 LET C4=.39 
75 LET 80=1.19 
80 LET T5=560 
90 LET M=1 

99 REM Q1=THICKNESS CORRECTION FACT0R 
240 PRINT ••T*',''L",”Y",'*P2","N” 

250 PRINT "V5","V6*’,"W8”,"W9",*’Y1“ 

253 PRINT •’L8*’,"H1 •*,*'Y2","W0" 

254 PRINT 

255 PRINT ••P5","P6",''’S'’,''P0","N0" 

257 PRINT '•D3",'*00","Fr',”F2*’,"F" 

258 PRINT '•Cr','*Q”,”Pl '•,"C0",’*K2" 

259 PRINT •’L0","L9" 

260 LET V=B 
270 LET D3 = ll 

280 PRINT "BELLOWS I.D.="D3;", DELTA V0L="V 

281 LET T=.010 
285 LET P2=25 

287 LET C0=2 

288 REM STRESS CORRECTED FOR THINNING 
290 LET F1=.445*F 

295 LET F2=F-F1 
299 LET Ql=l 


300 LET L=SQR(F1 *2*M*<T*Q1 ) T2/P2) 

301 IF L=>.15*D3 THEN 951 
310 LET D0=D3+2*L+2*CM*T-T) 

312 LET D2=D3+L+<M*T-T) 

315 LET A=<3.1416*D2t2)/4 

316 IF 01 <1 THEN 320 

317 LET 01 =D3/D0 

318 GO T0 300 

320 LET Y=V/A 

321 LET Z=.50 

322 LET Y1 =Z*Y 
325 LET Y2=Y-Y1 

387 REM Y5=T0TAL DEFL. F0R "Bl" V0iL 
328 LET Y5=B1/A 

330 LET N=3*Y2^E2*(T*01 )/CF2*2*.91 




T=600&750F 
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VANEP C0NTINUED 


350 LET W=2*3.14l6*D2*N*R*M*T»lc((L-n» @0„*T )_f3 . 1 41 6 6-0 1] 

360 LET Kl=3.1 416*E1 *D2*M*Tt3y(6*.a *N l^t3) 

365 LET K2=K1*E2/E1 

400 REM S2=S0L. HT. F0RME0 BELUbLjS 
440 LET S2=< <2*C0*T>+2*M+T3*N 
490 REM C1=FREE LENGTH 
500 LET C1=S2+Y1 

600 REM STRU; +.25' » XT FOP; jbTTlbCH. PP;gSS DOME: Z§*loTo 0F GttS DOME 

610 LET D8=D0+.25 

620 LET L8=(2*S2)+Y5+.80+2*.25 

624 REM OUTER GAN WEIGHT 

626 LET W8=3.l416=t=D8*.04*L8*R 

628 REM NAK DOME^ MOVE. HEAD, BELL. RINGS WEIGHT 
630 LET W8=W8+3.1416*R*(08t2/4*.06+CD0-.6)+.12*.6> 

632 LET W8 = W8+3.1416*R*C<D3-.26>t2/4*.06 + (D3-.26> + (S2+2 *-6o> *C)6> 

634 LET W8=W8+3.1416*R*D2*(L + .25)*.032>*=2 

636 REM W9=WT OF RESID NAK AT ZERO VOL DISPL AT lOOF 

638 LET W9=3.1416*(D2-.02)*L*(S2+.12)*54/1728 

640 LET W9=W9+3.1416*<D2-L>.10)+.10*<S2+.12)*S4/1728 

642 let W9=W9+D0t2/4*. 12*54/1728 

648 REM Z8=GAS DOME WEIGHT 

650 LET Z8=3.1416*R*(D0t2/4*.032+CD8*.032*.43 3 

652 LET Z8=Z8+3.1 41 6*D8*.08*.3*R 

660 LET W0=2*W+W8+W9+Z8 

680 REM V1=GAS VOL. IN ECU CAVITY 

685 LET VI =3.1416*CD2t2/4*(S2+Y5) + CD2-L-.32 t2/4* ( S3+. 62 3 5 

686 LET V5=V1 +3.1416*<C1/1 l*D0)t2*<15-<l/l 1 DO > /3 > -c 1 1 2*C 1 1 /3 3 3 : 

687 LET V5=V5+3.1 416*<D3-.06>t2/4*.7 

689 LET V5=V5+3.1 41 6*D0t2/4*.6 

690 LET V6=V5-A*Y 

691 LET P6=P2-(K2*Y2/A)*2 

692 LET P5=P6*V6*T5/(T6*V5) 

693 LET PI =P5-<K1*Y1/A)*2 

694 IF Pl=<4.05 THEN 697 

695 let P2=P2-.3 

696 GO TO 290 

697 IF PI =>4.0 THEN 700 

698 LET P2=P2+.05 

699 GO TO 290 

700 REM FAILURE MODE ANALi S=SECD VOL EXT TO BELLOWS 
705 REM PO=GAS PRESS, PRI BELLOWS FAILED, V0L=B 

708 REM NO=NAK PRESS, PRE BELLOWS FAILED, V0L=B 

710 LET S=3.1 41 6*<D2*L*.5*<Y5+2*S2)+( (D0+.06 )*.06*<Y5+2*S2+1 .75) ) > 

715 LET S=S+3.1416*D2*L*.88 

760 LET P0=T6*V5*P5/<T5*<V6+S+(V“B) ) ) 

775 LET N0=P0+2*K2*(CB-S)/A-Y1 )/A 
820 REM Q=SQUIRM PRESS AT FULL VOL 
825 LET 0=2*3.141 6*K2/(S2+Y5) 

830 REM H1=NATURAL FREQUENCY AT lOOF 
835 LET HI =9.85*SQRCK1 /W) 


« ft 


J 
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VtoNEP C0l\jriNUEO 


840 RSM L0=WELD LENGTH^ FEET 

845 LSI L0=C2*3.14l6*D0+.75+(L-2+C0*T+3.1416*C0*T)*2*N)/12 

846 LST L0-L0+<3.i416*D0)/12*3 

847 bET L9=L8+(1 .6/1 1 )^D0+.50 
900 SET DIGITS ^ 

910 PRINT T,L,Y#P2,N 
920 PRINT V5,V6, W8» W9#Y1 

930 PRINT L8>Hl>Y2>Wo 

931 PRINT 

935 PRINT P5,P6*S>P0,N0 
940 PRINT D3,D0.F1 .F2,F 
942 PRINT CUQ^Pl^Co>K2 

950 0)XINT Lq^L9 

951 PRINT 
999 END 


J 
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VANES 


2 REM SECONDARY M00P UNITS 

5 REM GAS BACK, |^EST I NG-Fo,RM 8ELLJ V=234 & 668, T=\oOQ*75oF 

6 REM 3 IDENTICAL LjNITBt 1 IN SEC, 2 IN PRI L00P 

7 REM BI=MAX CAPfl)3LE VOL 

8 LET 81=344 

9 REM 3«FAILURE M ^E V0m 

10 LET 8=236 
3q LSV T6=1q6o 
Sq let R=.297 

LET E1=28.8E6 
61? i^ST S2=26.61?S6 
70 LET C4=.39 
75 LET 80=1.19 
80 LET T5=560 
9o LET M»1 

99 REM Q1=THICKNESS CORRECTION FACTOR 
240 PRINT ••T","L","Y",’'P2","N” 

2$o PmSNT ”V5","V6'*, I(jo8* ,"W9","Y1” 

253 PRINT ••L8",”H1*%”Y8”,"Wq" 

254 PRINT 

255 PRINT "P5*',"P6'‘,"S*',"ft)0",”N0” 

257 JbRlNT ••03",'*D0", «F1 ■ ,”F2”,"F’' 

258 0)R3NT •'C1",”Q‘',*’P1*’,’ C0-',"K2” 

259 IDRcNT "L0","L9” 

260 LST V=8 
27q let D3=11 
275 LET T = =.qIO 

28q JDRINT "BELLOWS I.D.= 3l", DELTA V0L="V 

281 LET Y1 =1 .6153 

282 LET N=34.436 

283 LET L=. 52473 

285 LET V5=943.48 

286 LET P5=7.0502 
2t7 LST Cq=1? 

288 XSM STRESS CORlRSCTSD FOR THINlyjING 

289 LET &)0=D3+2 lfc+2*(M*T-T) 

290 LET 32=D3+L-|j(M+T-T) 

291 LSI A=<3.1416»D2r2>/4 

292 LET Y=V/A 

293 LET Y2=Y«Y1 
299 LET Ql=03/D0 

UQI IF L=>»15 &3 THSN 951 
3o2 i^ET V\o=V5“0) 

304 LET P6=P5*T6*V5/CV6*T5> 

327 REM Y5=T0TAL DEFL. FOR "81" VOL 

328 LET Y5=81/A 

350 LET W=2*3.1416*D2*N*R*M*T*(<L-2*C0*T>+3.1416+C0 ?3 
360 LET Kl=3.1416*El*D2*M*Tt3/<6*.91*C4+N*Lt3) 

365 LET K2=K1*E2/E1 

370 LET P2=P6+(K2*Y2/A>*2 
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VANES C0NTSNU^O 


375 
380 
385 
400 
440 
490 
500 
600 
610 
620 
624 
626 
628 
630 
632 
634 
636 
638 
640 
642 
648 
650 
652 
660 
680 
687 
693 
700 
705 
708 
710 
715 
7 60 
775 
820 
825 
830 
835 
840 

845 

846 

847 
900 
910 
920 

930 

931 
935 
940 
942 


LET F1=P2 ltt2/<2=)cMitCT*Ql)t2) 

LET F2=3 Yf-Z »8X'CT+0)1)/<2*.91 Bq 
LET F=F1+F2 

REM S2s=t0L. HT. FOITIMED BELLOWS. 

LET S2«<<2 90*T +2»M*T]*N 
REM C1=F><EE LEN TH 
LET C1=S2+Y1 

REM STRU: +.25" HT F0R ATTACH. PRESS D0ME: Z8=WT. 0F GAS DOME 
LET jb§=DQ+.25 

LET L8“<2 S2:+Y3+.80+2*.25 

REM OUTER CAN WEIGHT 

LET W8=3.1416*D8*.04*L8*R 

REM NAK D0ME# MOVE. HEAD, BELL. RINGS WEIGHT 
LET to8=Lo8+3.1416 W*<D8 t2/4^.06 + (D0- . 6 > ♦. 1 2*. 6) 

LET Lo 8 = WS+3. 141\o+X <K03-.26> t2/4*.06 + (D3-.26>*(S2+2*.60>*.06> 

LET W8=W8+3.1416 R &2 ^L+. 25)*. 032*2 

REM 0)9=:WT OF RSSIO NO)K AT ZER0 V0L DISPL AT lOOF 

LET lo9=3.1416 4D2 -.q2 3 *L* CS2 + . 1 2 >*54/ 1 728 

LET W9=W9+3.1416 «-D2-L- . 1 0 ) * . 1 0*CS2 + . 1 2 ) *54/1 728 

LET Io9-W9-HDOt 2/4 *12 gt-4/1728 

REM Z8=GAS DOME WEIGHT 

LET 78=3.1416 »<<30t2/4*.032+(D8*.032*.4) > 

LET Z8=Z8+3.1 16 »8 *q8*.3*R 

LET LoO &1 +Oi)8 + Ijo 9 +SJa 

REM V1=GAS VOL. IN ECU CAVITY 

REM GAS PRESS CALCS: P5/P6=INIT/FINAh <AT V0L=V3 Gibs' PR^SS 
LET Jbl =P5'-<K1 *Y1 /Jb)*2 C 

REM FAILURE MODE JbNtbL) S=SECD VOL EXT TO BELL0W_ 

REM P0=GAS PRESS, PRI BELLOWS FAILED, V8L=B 
REM N0=NJbK WRSSSi. JbRE BELLOWS FAILED, V0L=B 

LET S=3. 1416*(D2*L*.5*(Y5+2*S2)+( (D0+.06)*.06 4Yj+g fB+l.Tj])] 

LET S=S+3.1416*D2*L*.88 

LET P0=T6*V5*P5/CT5*(V6+S+(V-B>) > 

LET N0=P0+2*K2*< <B-S)/A-Y1 >/A 
REM Q=SQUIRM PRESS AT FULL VOL 
LET 0=2*3. 1416*K2/CS2+Y5) 

REM H1=NATURAL FREQUENCY AT lOOF 
LET HI =9.85*SQRCK1 /W> 

REM LO=WELD LENGTH, FEET 

LET L0=(2*3.1416*D0+.75+CL-2*C0*T+3.14l6 @q*V:*2 M3/12 

LET L0=L0+C3.1416*D0)/12*3 

LET L9=L8+(1 .6/11 )*D0+.50 

SET DIGITS 5 

PRINT T,L,Y,P2,N 

PRINT V5,V6, W8,W9, Y1 

PRINT L8,H1,Y2,W0 

PRINT 

PRINT P5,P6,S,P0,N0 
PRINT D3,D0,F1 ,F2,F 
PRINT Cl ,Q,Pl,C0,K2 


i 
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vo)N5!S continued 


950 PRINT l-0»L9 

951 PRINT 
999 END 
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